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Abstract 

We construct the complete set of metric-configuration solutions of 
the ghost-free massive bigravity for the scenario in which the o'—metric 
is the Friedmann-Lemaitre-Robertson-Walker (FLRW) one, and the 
interaction Lagrangian between the two metrics contributes an effec¬ 
tive ideal fluid energy-momentum tensor to the g-metric equations. 

This set corresponds to the exact background cosmological solution 
space of the theory. 

Keywords: Massive gravity, bigravity, self-accelerating cosmologies, 
dark energy. 
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1 Introduction 

The massive gravity theory which was constructed in BE] as a two-parameter 
family of actions, is a nonlinear generalization of the Fierz- Pauli massive 
gravity model [3]. In [2] it was shown that one of the actions of this theory 
is Boulware-Deser (BD)-ghost-free [H [5] up to the fourth order in metric 
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perturbations around flat space. An extension of this massive gravity model 
which originally admits a flat reference metric was also constructed in [HI El E] 
for a general background or reference metric. In [7], it was shown that the 
two-parameter family of actions given in |2] are BD-ghost-free at the complete 
nonlinear level for all perturbation orders. In |8], the analysis is extended 
to show that the general reference metric massive gravity theory in [B] is 
also BD-ghost-free at the complete nonlinear level for all orders. Later on, 
a ghost-free massive bimetric theory in which the interaction term between 
the foreground, and the backgronnd metrics arises from the mass terms was 
proposed by introdncing a copy of general relativity (GR) dynamics for the 
background metric [5]. Within this theory, a particular class of cosmological 
solutions of the conpled held eqnations of the two-metric-sectors PEHlEIlin] 
have been extensively studied in recent years [ElEslElESlEniETlllsllISl 
[20l 1211 |22l |23l [211 [25l [26]. The solution space corresponding to an effective 
decoupling of the two metric sectors of the theory is constructed in m 
The elements of the decoupling solution modnli that is derived in |27| give 
rise to self-accelerating cosmologies in the sector via the contribntion of 
an effective cosmological constant. On the other hand, in this work we will 
construct the general cosmological solution space of the massive bigravity 
theory. Thronghont our derivation, we will follow a parallel route with those 
of [28l [29] in which cosmological solutions of massive gravity are derived in 
a formal fashion withont explicit classification. On the contrary, in what 
follows likewise in [23) our main perspective will be to derive the entire 
set of metric conples (/, g) explicitly, which will consistently solve all the 
held equations of the theory, and which will enable a semi-deconpling of the 
/—metric from the gf—metric sector. Since we will consider a homogeneous, 
and an isotropic scenario in the sector apart from assigning a Friedmann- 
Lemaitre-Robertson-Walker (FLRW) form for g which corresponds to the 
cosmological background metric of the nniverse, we will also take the overall 
contribntion of the interaction Lagrangian to the metric eqnations as an 
effective ideal flnid energy-momentnm tensor. In this respect, the Bianchi 
identity of the interaction terms in the metric equations will natnrally 
transform to be the continnity eqnation of the effective flnid. We will show 
that when one proposes such a solution ansatz then the only remaining task to 
derive the set of metric couples which allow this picture is to And the general 
solntions of an inhomogeneous cubic matrix equation. We will derive the 
entire set of solutions of this matrix equation whose coefficients are fnnctions 
of the elementary symmetric polynomials of the solntions themselves rather 
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than being constants. By nsing the general solntions of this matrix eqnation 
one can constrnct the solntion moduli of the metrics (/, g) which admit a 
cosmological scenario in the sector. 

In Section one, starting from the bigravity dynamics by assuming a cos¬ 
mological sector and thus, introducing the above-mentioned ansatz for the 
contribution of the interaction terms in the metric equations we will de¬ 
rive an algebraic matrix equation whose solutions will lead us to the metric 
couples which are compatible with the held equations of the cosmological 
picture. In Section two, we will derive the Jordan normal form solutions of 
this equation. By using this complete set of Jordan normal form solutions, 
and other three sets constructed from them we will discuss in Section three 
that when special form of similarity transformations are used one can obtain 
the general set of solutions of the above-mentioned ansatz-generated matrix 
equation. An outline of the proof of the completeness of this set of solutions 
will be given in the Appendix. Then, again in Section three as a consequence 
of the completely-derived general solution space of the ansatz equation we 
will be able to dehne the cosmological solution moduli space of bigravity. 
We reserve Section four for the derivation of the equations of the cosmo¬ 
logical dynamics in sector, as well as a discussion about the associated 
/—dynamics, and its solution methodology. We will also present the outline 
of the /-sector solution construction of an example. 


2 The set-up 

The ghost-free bigravity action can be given as 


^ = - 


IbvrG 

K 

IbvrG 


J 


R3 + 


+ Si 


M 


dx^ 




+ 


dZTD 


where g is the foreground, and / is the background metric which are coupled 
to two types of matter via the actions respectively. A^, A-f are the 

cosmological constants in each sector. R^, R^ are the corresponding Ricci 
scalars. The Lagrangian which describes the interaction between the two 
metrics above reads 


= /3iei(\/S) + /? 2 e 2 (\/S) + /? 3 e 3 (\/S), (2.2) 
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where {e„} are the elementary symmetric polynomials 

Cl = ei(\/S) = 

62 = e2(\/S) = ^((tr\/S)^ - tr{'/Ef), 

63 = e 3 (\/S) = ((tr\/S)^ - 3 tr\/S tr(\/S)^ + 2 tr(\/S)^), (|2.3p 

0 

corresponding to the square-root-matrix 

= ypu. (2.4) 

Originally, the interaction Lagrangian also contains the terms PoCq = /3q, 
and 13464 = /94det\/S. However, we combine their contributions with the 
cosmological constants A^, and Af, respectively. Eq. fl2.2p reduces to the 
Fierz-Pauli form in the weak-field limit when one chooses |H] 

f3i + 2/?2 + (33 = —1. (2-5) 

From Eq. fl2.ip one can obtain the held equations for the metric g as 

( 2 . 6 ) 

Whereas, the held equations of the metric / become 

'‘[A. - + ’"At = (2.7) 

The corresponding energy-momentum tensors arising from the interaction 
term Eq. fl2.2p are 

= -9^lp'r% + ( 2 . 8 ) 

and 

pi = (2.9) 

The matrix r with the entries {r^^} is dehned to be na 

r = /? 3 (\/E)^ — {(32 + /? 3 ei)(\/E)^ -|- (/?i -1- /? 2 ei + I3z^2)^^- (2-10) 

The ehective energy-momentum tensors in Eqs. fl2.6p . and fl2.7p are ought to 
satisfy the Bianchi identities 

= 0, v;(r 0 ". = 0. (2.11) 
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If a solution configuration satisfies one of these equations then the other one 
is automatically satished [T5[ [TB] . Now, let us focus on the cosmological 
solutions in the sector of the action Eq. fl2.1l) . Thus, we will take g as the 
FLRW metric 

g = —dt^ H— + a^{t)r‘^d6‘^ + a^{t)r‘^sin^6dip‘^. (2.12) 

1 — kr^ 

Let us also consider the solutions for which the effective energy-momentum 
tensor entering into the metric equations in Eq. fl2.6p takes the form 

= (Kt) + P{t))U^U^ + p{t)gf,^, (2.13) 

of an ideal fluid. For an ideal fluid the on-shell Lagrangian can be taken as 


= P, 


(2.14) 


so that 


(2.15) 


Therefore, for the solutions which generate Eq. fl2.13p in the held equations, 
Eq. (12. 6 p a close inspection of the interaction term in the action Eq. (12.ip . 
and Eq. (I2.15p shows us that we must have 


k^int P- 


(2.16) 


This is due to the fact that while Eq. fl2.8p is derived by varying the inter¬ 
action term in the action Eq. fl2.ip with respect to the metric g, Eq. fl2.15p 
is obtained by varying the Lagrangian Eq. fl2.14p with respect to g, and by 
using the hrst law of thermodynamics [29]. Our effective huid that is intro¬ 
duced in Eq. fl2.13p will certainly obey the hrst law of thermodynamics as, 
it must satisfy the conservation equation in Eq. fl2.1ip which will result in 
an ordinary continuity or huid equation when Eq. fl2.13p is substituted in it. 
If we take the ehective ideal huid four-velocity vector as = (1,0, 0,0) in 
the rest frame of the huid, and use the FLRW metric Eq. fl2.12l) we obtain 


g-^T^ 


j-p 0 0 o\ 
0 p 0 0 
0 0 p 0 

y 0 0 0 pj 


(2.17) 
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Index raising on both sides of Eq. fl2.8p by the metric g gives 


(T^) M 

(2.18) 

where i'l h = [a~^'l ^Jh.. Bv usine; Eqs. 02.1(111. and 02.17P 
we obtain the matrix equation 

in this expression 

+ B{VEf + C'\/S + V = 0, 

(2.19) 

with 

I'D 0 0 0\ 

^ 0 0 0 0 

0 0 0 0’ 

\0 0 0 0/ 

where 

(2.20) 

B = —^2 — 

C = fXl-X P 2^1 + (33^2, 

D = —p — p. 


3 Classification of the solutions 


Next, we will derive and classify the Jordan canonical form solutions of the 
cubic matrix conation (|2.19p. This is a hiehlv non-trivial matrix conation for 
two reasons: hrst, it is not in a polynomial form, and second, its coefficients 
are functions of the elementary symmetric polynomials Ci, 62 of its solutions 
\/S rather than being constants. For this reason, in this section we will derive 
the diagonal and the nondiagonal Jordan form solutions of it, and show in 
the next section that they can be used to generate the entire solution space. 
Firstly, let us define the polynomials 

Ax^ + Bx^ + Cx + D = 0, 

(3.1) 

and 

x{Ax^ + Bx + (7) = 0, 

(3.2) 


whose roots we will generally call a*, and Xj, respectively. Since, in general 
for any Jordan canonical form matrix J (diagonal or nondiagonal) when it is 


6 













substituted into the Eq. 02.191) the eigenvalues namely the diagonal elements 
of J must satisfy one copy of Eq. 03.ip . and three copies of the polynomial 
in Eq. 03.2p the multiplicity of a* in the diagonal of J must be one. This 
is obvious for the diagonal Jordan forms. Besides, if we have a nondiagonal 
Jordan canonical form as 


(ai 1 0 

0 ai • 


\0 0 0 

when this solution ansatz is used in Eq. 02.19p then the diagonal entries 
would lead to two inconsistent ecmations A{aiY + + C{ai) + D = 0, 

and A{aiY + + C{ai) = Cu. Therefore, following this observation we 

can conclude that the Jordan form solutions of Eq. fl2.19p are partitioned as 


0 \ 


(O 


J 


Oii 

0 

0 

i^3x3 


(3.4) 


where iJsxs is a three by three Jordan normal form matrix which satishes 
the matrix polynomial equation 


+ C'hfsxs — 0. (3-5) 


The diagonal ones are naturally in the form of Eq. fl3.4|) . and the nondiagonal 
ones must be in this form due to the partitioning nogo fact discussed above. 
7^3x3 must be constructed by its eigenvalues {A^}. As a result, we deduce 
that the classihcation of the Jordan form solutions of Eq. 02.191) must be 
based on the classihcation of the solutions of Eq. 03.Sp . We should remark 
one important point here that, in deriving these solutions we will exclude 
the cases which arise from the conditions B = 0, and/or (7 = 0. As it will 
be clear in Section four, such conditions would lead to extra constraints on 
the equation of state p = p{p) of the effective huid that would cause it to be 
nondynamical. In addition, we will also exclude the trivial case of = 0. 

^Obviously, these two equations can be satisfied simultaneously ii D — 0 which would 
correspond to r = 0 cases. However, we take D 0 so that the effective fluid is not simply 
an effective cosmological constant. Thus, we exclude the r = 0 solutions here as they are 
completely derived in m- 
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3.1 A > 0 Solutions 


Let us first consider the solutions when A = B'^ — AAC > 0. In this case there 
are three distinct roots of the polynomial in Eq. f|3.2p . they are {0, Ai, A 2 }. 
Now, none of the Jordan normal forms which satisfy Eq. fl3.5|) may have a 
repeated root of its minimum polynomial which is formed by a subset of the 
factors in Eq. fl3.5l) . That is to say, for these cases if we write Eq. fl3.5p in 
the form 

H{H = (3.6) 

where I 3 is the unit 3x3 matrix then we see that the solutions of this equation 
must make the product of three factors, or any two factors, or just a single 
factor vanish. Thus, there are 1 + 3 + 3 distinct classes of solutions where 
one of them is trivial. For each class the vanishing combination of factors 
become the minimum polynomial of the corresponding 3x3 matrix solution. 
Therefore, we also observe that for each of these classes the roots of the 
corresponding minimum polynomial are not repeated (as Ai, A 2 are distinct ) 11 . 
Since a matrix is diagonalizable if and only if its minimum polynomial has 
no repeated roots we deduce that in this case any solution satisfying Eq. 
fl3.5l) must be a diagonalizable one. On the other hand, if a Jordan canonical 
form satishes Eq. fl3.5p its similarity equivalence class also does. Hence, 
combining these two facts we conclude that all the Jordan forms which satisfy 
Eq. fl3.5p must be diagonal. The reader may also verify this result by direct 
substitution. In other words, none of the nondiagonal Jordan canonical forms 
whose diagonal elements are chosen from the set {0, Ai, A 2 } satishes Eq. fl3.5p . 
Therefore, upon this identihcation of the 3x3 sectors, the corresponding 4x4 
Jordan normal forms which satisfy Eq. fl2.19p can be listed as 


(ai 

0 

0 

0\ 


(at 

0 

0 

0\ 


/ at 

0 

0 

0\ 

0 

Ai 

0 

0 

,iV2 = 

0 

A 2 

0 

0 

,iV3 = 

0 

0 

0 

0 

0 

0 

Ai 

0 

0 

0 

A 2 

0 

0 

0 

Ai 

0 

^0 

0 

0 

aJ 


VO 

0 

0 

A2y 


VO 

0 

0 

aJ 


^We should note that, if Ai, or A 2 vanish then there will be a repeated root (the zero 
root). In this case a single nondiagonalizable Jordan form exits but this case requires 
C = 0 condition that we exclude as we pointed above. 





















^ 0L{ 

0 

0 

0\ 



0 

0 

0\ 


^ Oil 

0 

0 

0\ 

0 

0 

0 

0 

, iV5 = 

0 

0 

0 

0 

, N,= 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A 2 

0 

0 

0 

0 

0 

Vo 

0 

0 



VO 

0 

0 

X 2 } 


Vo 

0 

0 

A 2 / 


(IS2I) 


as well as the ones, 



^ Oil 

0 

0 

0\ 



^ Oil 

0 

0 

0\ 

w = 

0 

Ai 

0 

0 

,iV8 


0 

A2 

0 

0 

0 

0 

Ai 

0 

— 

0 

0 

A2 

0 


Vo 

0 

0 

A 2 y 



Vo 

0 

0 

Ai/ 





^ Oil 

0 

0 

0\ 






IV9 


0 

Ai 

0 

0 






— 

0 

0 

A2 

0 








Vo 

0 

0 

V 





(M 


in which both of the roots Ai, A 2 appear. By direct snbstitution, the reader 
may verify that these matrices do solve Eq. fl2.19p . and the corresponding 
nondiagonal Jordan forms that share the same eigenvalues can not satisfy 
Eq. fl2.19p in this case when A > 0. To hnd the explicit form of these 
solutions we have to know ei, 62 which constitute both the coefficients given 
in Eq. fl2.2ip (of the matrix equation that these solutions must satisfy), and 
the entries of these solution matrices listed above. In other words, we have 
to solve ei, and 62 in terms of the {/dj}— parameters of the action Eq. fl2.ip . 
and the constituents of the solution ansatz Eq. fl2.13p so that Eq. fl2.19p is 
satished. We will hrst consider the cases iVi^ 2 , 3 , 4 , 5 , 6 - If we take the trace of 
these solutions we get 


ei = tr{Na) = ai + nXj, (3.9) 

where n = 3,3, 2, 1, 2, 1 for Ni, N 2 , N^, N 4 , N^, Nq, respectively. We also have 
tr(W)2 = (ei)2 - 262 = («.)' + (3.10) 

By using Eq. fl3.9p . and singling out 62 from this expression we get 


62 = 


n + n 


(Aj)^ + neiXj. 


(3.11) 
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(3.12) 


If we substitute this result into Eq. fl3.2p we obtain the relation 

fl(Aj) “1“ c = 0, 


where 

f2 — Tl — 

a = - - -/^s, b = (n — l)/? 3 ei — (32, c = (3i + /? 2 ei- (3.13) 

Since from Eq. fl3.9p we have a* = ei — nXj, substituting this into Eq. fl3.ip . 
and successive usage of Eq. fl3.12p leads us to the relation 


A, = 


(n^ - u)/?3 


b 

-ei 

a 


b c\ ^ o 2b o 

n— +n- \ + P 2 nei + (32n - (3in 

a / a 


-1 


X 


(n^ — n)(3^cei ^ n(n^ 


n)(33bc rPc(32 


2 a 


2a2 


- (3iei - D 


(EH 


Finally, when we use Eqs. fl3.14p . and fl3.1ip back in Eq. fl3.12p and we refer 
to the dehnitions in Eq. fl3.13p we obtain an equation for ei solely in terms 
of the /9j—coefficients. For the n = 3 cases this equation reads 


®3(ei)^ + ^ 3 ( 61 )^ + C 3 (ei)^ + ^361 + /3 — 0, (3.15) 


where we dehne 


03 = 3(/33)^(—3(/?2)^ + d/^i/^s), 

^3 = 6/33(—15(/32)^ + 20/3i/ 32/33 + 2(/?3)^iA), 

C 3 = —216(/32)^ + 159/3i(/ 32)^/33 + 172(/3i)^(/33)^ + 102/32(/33)^Zl, 
ds = 204/3i/32(—3(/?2)^ + 4:/3il3^) + /33(249(/32)^ + 20/3i/33)iA, 

/a = -432 (/?i)2(/32)' + 57Q{^if(3s + 3Qi^2fD + 2A0(3MD - 125{(3sfD\ 

(EH 

For the n = 2 cases we get 

^ 2 (^ 1 )^ + ^ 2 ( 61 )^ + C 2 ei d2 = 0, (3.17) 

where 

O 2 = —(/^2)^/^3 + ldi(32{.(3zY 1 

&2 = —6(/?2 )‘^ + 4/3i(/32)^/53 + ‘^{l3iY{l3zY + I32{I3zYD, 

C2 = -2iPi{P2f + 21(A)'/92/93 + KP2fP^D - 2Pi{P^fD, 

d2 = —73{(3iY{(32Y + 18(/^i)^/^3 + 3(/32)^-D + Q(3i(32(3sD — 4(/33)^iA^. (|3.18p 
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For each real root of Eq. fl3.15p we have the solutions Ni,N 2 , and for each 
real root of Eq. 03.171) we have the solutions of Eq. 02.19P with the 

corresponding entries that can be read from 


A,= 


-b ± — 4ac 

^ ’ 


a* = Cl — nXj 


(3.19) 


The domain of validity of these solutions are determined by the conditions 
— AAC > 0 , and 6 ^ — 4ac > 0 , together with the reality conditions of the 
corresponding roots of the Eqs. 03.15p . and 03.17P which can be obtained 
from the dehnitions of the coefficients in Eqs. fl3.16p . and fl3.18p . These 
conditions will dehne a validity domain for each particular solution in the 
union of the parameter space, and the state space of the effective ideal 

fluid. On the other hand, for the cases with n = 1, namely for the solutions 
of the form Nq we have a simpler picture. In these cases, Eq. fl3.12p gives 


Xj — Cl + 


A 

/92’ 


(3.20) 


and from Eq. fl3.1ip we have 


62 = -(ei + + ei(ei + 

P2 

Now, by using Eq. in Eq. ([33]) we get 


tti = 


A 

132 ' 


Substituting this result, together with Eq. (I3.2ip 

/52 



(3.21) 


(3.22) 

into Eq. (13. ip gives us 

(3.23) 


The Eqs. fl3.20p . fl3.22p . and fl3.23p dehne the explicit form of the entries of 
N 4 , Nq in terms of the parameters of the theory, and the state of the effective 
huid. The validity domain of these solutions is dehned from the condition 
B^ — 4AC > 0 without extra requirements. Now, let us consider the solutions 
A^ 7 , Nq. Taking the trace of these solutions we get 


ei = trNa = ai + 2Xj + Xk 


, B 

ai + Xj - —, 


(3.24) 


11 




































where we label the excess or the repeated root on the diagonal of the solution 
by Xj. By refereing to the definitions in Eq. fl 2 . 2 ip we find that 

= -A, - (3.25) 

/?3 


We also have 


tr{N,y = {eif - 262 = + 2(A,)2 + (Afc)l 

By using the identity 

(Ajj + [Xk) - ^ - 

in the above equation we see that for these solutions 

„ , _ A 
‘'“W A' 

By using this result in Eq. fl3.25p we obtain the equation 


/A n2 /^2 , Pi {P2Y 


= 0 , 


for Ai. Its solutions are 


± ^ (P 2 ^ 


.(ML 




(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 


Using Eqs. (I3.29p . and (I3.3np in Eq. fl3.2p will enable us to write 62 in terms 
of 61 . After some algebra we get 


62 = 


2/93 2 


_3(M + 4h 

(hf h 


eiT 


Ps 





(3.31) 


We note that, when the (+) solution is taken in Eq. fl3.30p the opposite 
sign must be chosen in Eq. fl3.3ip . and vice versa. Finally, substituting 
expressions fl3.25p . fl3.30p . and fl3.3ip into Eq. fl3.ip will lead us to the explicit 
value of 61 for either of the solutions in Eq. (I3.3np . For the solutions A^ this 
computation reads 


± ^ _ 2{PspD 

' P3 (U± 


(3.32) 
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where we defined 


= 3(/?2)% - ± m^?\ 


In summary, for these latest cases we find that 








(3.33) 



(»t 

0 

0 

0 \ 



0 0 

0 \ 


0 

A+ 

0 

0 


0 

A” 0 

0 

Nj — 

0 

3 

0 


0 

1 Ng — 

0 

J 

0 AT 

0 


VO 

0 

0 

k) 



0 0 

\ / 

are the solutions 

of Eq. 

(EH- 

Explicitly, together with 

Eqs. 


(E31 


(13.32^ we have 


af = -A± 


^ + — 


A± 


(3.35) 


The domain of validity of these solutions in the parameter space of the action 
Eq. (12.ip . and the state space of the ansatz Eq. (I2.13p is governed by the 
conditions B'^-AAC > 0 , and - 3 (/ 92 )^/(/d 3 )^+ 4 /Si //53 > 0 with the respective 
substitutions of 62 , and Ci from Eqs. (13.311) . and (13.321) . The final solution 
we have to derive explicitly in this class is N^. If we take its trace we find 
that 

€-\ = trNg = Oj T Ai T A2 = ttj — 


A 

From this relation by referring to Eq. (I2.2ip we see that 

Also, 

triNgf = {e,Y - 262 = («.)^ + (Ai)^ + (A2)^ 

By using the identity (I3.27P this relation reduces to the condition 

(/92)'-/5i/33 = 0. 


(3.36) 

(3.37) 

(3.38) 

(3.39) 


Now, substitution of Eq. (I3.37p into the polynomial Eq. (13.Ih gives 62 in 
terms of ei explicitly. It reads 


62 — 61 

P3 


/di D 


(3.40) 
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We see in this formulation that e\ remains completely an arbitrary spacetime 
held. For a particular choice of it one can read from Eq. fl3.4Up . and 
construct A,B,C explicitly in terms of {A}, and D via their dehnitions in 
Eq. fl2.2ip . then one can explicitly obtain the entries of Ng from 


A 


1,2 — 


-B ± - AAC 

2A 


(3.41) 


and Eq. (I3.37p . For this solution to exist the conditions Eq. (I3.39p . and 
B^ — 4AC > 0 must be satished. Again, the second of these dehnes a domain 
in the union of the action-parameter space of the theory, and the state space 
of the effective huid. 


3.2 A = 0 Solutions 

We now turn our attention on the cases when A = B^ — 4AC = 0. In these 
cases there is a repeated root A' = —B/2A of the polynomial Eq. fl3.2p . 
The roots of Eq. (13.2p become {0, A', A'}. Since, when it is factorized Eq. 
(13. 5 p has a repeated factor some of the Jordan normal forms which satisfy 
Eq. (13. 5 p may have a repeated root of their minimum polynomials. Thus, 
when A = 0 we have nondiagonal, as well as diagonal Jordan normal forms 
which satisfy Eq. fl3.5p . We again, do not consider the solutions arising from 
i? = 0, and/or C = 0 conditions. Bearing this fact in mind, therefore, in this 
case the list of all the possible Jordan normal forms which satisfy Eq. (I2.19p 
can be given as 


(oii 

0 

0 

o\ 


(Oli 

0 

0 

0\ 


/ Oii 

0 

0 

0\ 

0 

A' 

0 

0 

,iVii = 

0 

A' 

0 

0 

,iVl2 = 

0 

A' 

0 

0 

0 

0 

A' 

0 

0 

0 

0 

0 

0 

0 

A' 

0 

VO 

0 

0 

^7 


VO 

0 

0 

V 


Vo 

0 

0 

0/ 


/ 

0 

0 

0\ 


^ Oii 

0 

0 

o\ 

0 

A' 

1 

0 

, iVl4 = 

0 

A' 

1 

0 

0 

0 

A' 

0 

0 

0 

A' 

0 

VO 

0 

0 

^7 


Vo 

0 

0 

0/ 


(E321) 
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The reader may again verify that these matrices satisfy Eq. fl2.19p by direct 
substitution. If we take the trace of the matrices in Eq. fl3.42p we find 


Cl = trNa = tti — n 


B 

2A'' 


(3.43) 


where n = 3,1,2, 3, 2, for A^io, iVn, iVi 2 , iVi 3 , A^i 4 , respectively. When the 
definitions in Eq. fl2.2ip are used this relation yields 


n (32 2 — n 

Oii — -;;—h -Cl. 

2/^3 2 ' 


We also have 


B 


tr{N,f = {e,f-2e2 = {a.f + n{—)\ 


which can be written in the form 

n^ + n(/l 2 )^ in? — n (32 n? — ?>n . ,3 


62 = 


8 


/Ss 


ei - 




(3.44) 

(3.45) 

(3.46) 


by using Eqs. fl2.2ip . and fl3.44p . Now, substituting Eqs. fl3.44p . and fl3.46p 
into Eq. fl3.ip yields 


a'(ei)^ + &^(ei)^ + c'ci + d' = 0, (3.47) 


where 


a' = (— 2 n + 3n^ — in?){(3-if, 

h' = (4n + (32{(3zY, 

d = (16 - 9!n)(3i{(3df + ( 6 n - 3n^ - 3n^){(32^(33, 

d! = -{3rd + rd){(32Y ~ (|3.48|) 

On the other hand, for these solutions to exist we have the condition A = 
B'^ — 4AC = 0. Again, by using the definitions in Eq. (I2.2ip . and by 
expressing 62 in terms of ei via A = 0 condition, then by substituting the 
result in the equation fl3.46p we obtain 


n 


3n + 2 


8 


(ei)^ + 


n 


-n-2/32 n‘^ + n + 2(/32y 




ei + 


{dzd 




(3.49) 


For n = 1, this equation is reduced to a linear one and it has the solution 


ei 


Ml , ^ 

d2 dd 


(3.50) 
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for n = 2, the ci—terms vanish and it boils down to the condition 


{^ 2 ? - /Si/Ss = 0, 

and for n = 3, it has the solutions 

2ft /ft 3 (ft)2 

ft V* "“W 


(3.51) 


(3.52) 


provided that /^iZ/Ss — 3(/92)^/4(/93)^ > 0. When Eqs. fl3.50p . and fl3.52p are 
substituted into the equation fl3.47p one hnds the equation of state for the 
effective fluid. However, at this stage we need not explicitly derive these 
expressions since such solutions can not lead to evolving scale factors for 
these cases. We will explain the reason why this occurs for the n = 1, and 
n = 3 cases, and thus, why we can disregard them in the next section. On 
the other hand, for the n = 2 cases namely for the solutions W 2 , -^14 Eq. 
fl3.49p does not £x ei in terms of the /d*—coefficients but only results in a 
condition on them. For these cases Ci must be solved from Eq. fl3.47p thus, 
the equation of state of the effective ffuid is not hxed. Solving Eq. (I3.47p 
which reduces to be a quadratic equation when n = 2 for ei yields 


e 


± 

1 


—c' ± 


26' 


n=2 


(3.53) 


Having found ei now, we can explicitly express the entries of the solutions 
Ni 2 , Ni4 via 




A' = I 

^ 2 A 2/^3 2 ’ 


= “TT 


/Ss’ 


(3.54) 


where Ci must be substituted from Eq. 03.531) . We see that, there are two 
sets of solutions for each of iVi 2 , and Nu. For the existence of these solutions, 
there are two conditions to be satished; one of them is Eq. 03.5ip . and the 
other one is (c')^ — 46'd' > 0. 


3.3 A < 0 solutions 

When A = —4AC < 0 the polynomial fl3.2p has complex roots. The roots 
of Eq. fl3.2p are {0, A, A*}. Upon factorization, Eq. fl3.5p has complex root 
factors too. The non-trivial Jordan normal forms which satisfy Eq. fl3.5l) by 
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causing a minimum polynomial that is a sub-factor in the factorization of Eq. 
fl3.5p to vanish must have complex eigenvalues. Thus, they are nondiagonal. 
In this case, if A is an eigenvalue (or a root of the corresponding minimum 
polynomial) of any Jordan normal form which satisfy Eq. fl3.5p A* must also 
be an eigenvalue. Beside this fact, by assuming B ^ 0, and C ^ 0, also by 
considering the form in Eq. fl3.4p we conclude that, the only possible Jordan 
normal form in this class that would satisfy Eq. fl2.19p is 


A^15 = 


/a* 

0 
0 
VO 


0 0 o\ 

RIO 
-IRQ 
0 0 0 / 


(3.55) 


where we dehne \ = R + li with 

—4 


/ = 


^JAAC - 52 
2A ' 


(3.56) 


Again, it can be verihed that Eq. (I3.55p satishes Eq. (I2.19p via direct 
substitution. Now, by taking the trace of A^is, and also by referring to Eq. 
fl 2 . 2 ip we hnd that 

a,; = (3.57) 




Furthermore, 


= (ei)2 - 262 = (a.)' + 2{R^ - J^), (3.58) 


which reduces to the condition 


{^2? - m = 0, 


(3.59) 


upon using the dehnitions in Eq. 
fl3.1ll leads us to 

2/32 


62 = 


-61 - 


fl2.2ip . Substituting Eq. fl3.57p into Eq. 

m? 


A 

/Ss 


- vA + 


D 


(3.60) 


(/Ss)^ /Ss ' /32' 

We realize that, in this solution 61 remains to be an arbitrary spacetime 
function. When one specihes 61 one can express 62 in terms of it from Eq. 
fl3.60p . then one can explicitly obtain the matrix entries of Ws by using Eq. 
fl2.2ip in Eq. fl3.56p . and from Eq. fl3.57p . The conditions of existence of iVi 5 
are Eq. fl3.59p . and A = — 4AC < 0 in which the particular choice of 61 , 

and the corresponding 62 must be used. 
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4 The solution space 

In the previous section, we have explicitly constructed the entire set of 
nontrivial Jordan canonical form solutions of Eq. fl2.19p . We have disre¬ 
garded the trivial case of \/S = diag(Q!i, 0 , 0 , 0 ) which results in nonphysical 
/—metric solutions. Before dehning the solution space of Eq. (I2.19p . let us 
discuss one last constraint on the solutions that we constructed in the previ¬ 
ous section. In obtaining the Jordan normal form of the solutions, although 
we used the conditions on the elementary symmetric polynomials ei, and 
we did not refer to the 63 —structure of the solutions. This is a necessary, and 
a crucial point, as our solution ansatz Eq. fl2.13p brings the constraint Eq. 
fl2.16p on Cl, 62 , 63 values of the solutions which we have not yet considered. 
To impose this condition on the solutions we simply take the trace of Eq. 
fl2.19p . Following the trace operation on Eq. fl2.19p if we make use of the 
Eqs. fl2.3p . and (12.16p we get the relation 


— 2/?iei — ^• 2.^2 — p J- 2p — 0. (4.1) 

When, for each solution the appropriate value of e\ = ei{(3i,p,p), and 62 = 
e 2 (A) PyP) of fhcit particular solution are used in Eq. (14. ip the above expres¬ 
sion fixes the equation of state of the effective fluid that is p = p{p) for the 
solution chosen. At this stage, we can explain why we have to exclude the 
solutions of Eq. (I2.19P for the B = 0, and/or C = 0 cases. These conditions 
via the dehnitions in Eq (I2.2ip will bring an extra constraint on p, p. Thus, 
when solved simultaneously with Eq. (14. ip this constraint will cause p, and 
p to be constants. Besides, the n = 1, 3 cases of the A = 0 solutions of the 
previous section also can be eliminated as cosmological solutions. Similarly, 
for those solutions we have seen that the A = 0 condition already caused the 
determination of the equation of state of the effective fluid composing the 
cosmological solution ansatz Eq. (I2.13p . For these cases, when one solves the 
resulting conditions on p,p coming from Eqs. (I3.47p . and (l3.5Up . or (I3.52p 
together with Eq. (14. ip one sees that both p, and p must be constants again. 
Therefore, as it will be clear in the next section for all of these cases the scale 
factor can not evolve hence, it results in static, and nonphysical cosmological 
solutions. Now, hrstly let us dehne the set of Jordan normal form solutions 
ofEq. (ICTD 

Jl = |iVi,iV2,iV3,iV4,iV5,iV6, iV7,iV8,iV9,iVi2,iVi4,iVi5|. (4.2) 
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Next, we introduce the matrix field 


) 

P3x3{xn 


(4.3) 


Pl = 


m{x^) 


where P{x^) is an invertible 3x3 matrix field, and m{x^) is a scalar field 
which can simply be taken as m(x^) = 1 without loss of generality. Since 
any element J G J7i is a solution of Eq. fl2.19p if we perform a similarity 
transformation on both sides of Eq. fl2.19p we get 

A{P^^JPi)^ + 5(Pf VPi) 2 + CPf VPi + P = 0, (4.4) 

where we have used the fact that P^^PPi = P. Therefore, we see that 
Ph^ JPi is also a solution of Eq. fl2.19p for any J G jTi, and for any matrix 
field of the form Eq. fl4.3p . In that regard, we can define a subset of the 
solution space of Eq. fl2.19p as 


Ml = 


|\/S I \/E = PfVPi 


VJ G J7i, and detPi ^ 0 


(4.5) 


in which Pi is any matrix field of the form Eq. fl4.3p . We should remark that 
Eq. fl4.ip that is obtained for the element J G J7i remains the same for also 
the corresponding element Pi“^ JPi G Mi as the elementary symmetric poly¬ 
nomials do not vary under similarity transformations. Now, let us consider 
the matrix equations 


A{jf + B{jf + CJ + Vi = Q, (4.6) 

where i = 2,3,4. Here, P 2 = diag(0, P, 0, 0), P 3 = diag(0, 0, P, 0), and 
P 4 = diag(0, 0, 0, P). The Jordan canonical form solution spaces namely JT) 
of these equations can be constructed from the elements of Ji. In particular, 
for example, the diagonal elements of J 2 are obtained by placing to the 
second diagonal entry, and by shifting the rest of the diagonal entries diag¬ 
onally in the elements of jTi. Also, the two nondiagonal elements of are 
obtained again, by placing a, in the second diagonal entry and by shifting 
the primary blocks diagonally in Nn^Ni^. The elements of jTs, and J 4 can 
be obtained in a similar fashion. We should state that all these diagonal 
shifting operations which are used to generate the elements of 1 / 2 , 3 , 4 from the 
elements of J7i do not change the elementary symmetric polynomials of the 
corresponding element as we keep the diagonal content in these operations. 
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Thus, the parametrization derived in Section three for the elements of J7i are 
also valid for the elements of 1 / 2 , 3 , 4 . Next, let us define the invertible 4x4 
transformation matrix functions 


(0 

^2 

0 

0\ 


/o 

0 

ts 

0\ 


(0 

0 

0 

u\ 

• 

0 

• 

• 

,^3 = 

• 

• 

0 

• 

,^4 = 

• 

• 

• 

0 

• 

0 

• 

• 

• 

• 

0 

• 

• 

• 

• 

0 

V* 

0 

• 



V* 

• 

0 



V* 

• 

• 

0/ 


where ^ 2 , 3,4 are arbitrary functions of like m{x^), and the entries which are 
not specified in the above 4x4 matrices form up partitioned 3x3 invertible 
matrix functions. If now, we define the solution spaces 


A42= j 

[^1 

x/s 

= P2JP2" 1 

VJg 

J2, 

and 

detP2 

# 0 ] 


-M3= ] 

[^1 

x/s 

= PzJP^^ 1 

VJ e 

iTs, 

and 

detP3 

^ 0 ] 


A44= j 

[VEI 

^/s 

= PaJP^^ I 

VJe 

Xa, 

and 

detP4 

^ 0 ] 











(m 



M -- 

= A4 1 U Xi"} 

,uA43 

UM 4 , 



(4.9) 


becomes the general solution space of the Eq. (I2.19p . We will give a sketch of 
the proof of this fact in the Appendix. Now that we have found the complete 
solution space of Eq. fl2.19p . we can turn our attention on the background 
metric solutions of the action Eq. fl2.ip . By referring to Eq. fl2.4p we now 
have 

/ = (4.10) 

where S is the square of any element in A4, and g is the FLRW metric. 
However, not all elements of M. which solve Eq. fl2.19p will lead to symmetric 
results in Eq. fl4.10p thus, physically acceptable background metrics. We 
have to impose the condition 

gE = E^g, (4.11) 

which guarantees the symmetry of /. Therefore, we define the cosmological 
solution moduli of the action Eq. fl2.ip as the set 

rc = {{ 9 j) \ f = 9 X^ \XeM,e.ndgX^ = {X^yg}. (4.12) 
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In special, when one chooses the diagonal elements in J 2 -, Jz-, then 
squares them, and substitutes the result in Eq. fld.lUp one directly obtains 
the exact background metric solutions in a concise way without being obliged 
to concern the symmetry condition. On the other hand, for the more general 
cases one has to choose a special form for the matrices Pi, P 2 , -P 3 , Pa to satisfy 
Eq. fld.lip . Since, the symmetry requirement in Eq. fl4.1ip becomes 

gPiJ^P-^ = {P,J^P-Yg, (4.13) 

a closer inspection denotes that for a particular choice of J G JT) this equa¬ 
tion brings three algebraic constraint conditions on the function-entries of 
the solution-generating P—matrices (in particular, their unspecihed 3x3 
partitions) which enable us to determine three of the entries of these parti¬ 
tions in terms of the other six entries which remain arbitrary. Next, we will 
give a summary of the cosmological dynamics. 


5 Cosmological dynamics 


In the sector beside the effective ideal fluid energy-momentum tensor that 
is introduced in Eq. fl2.13p . we will also take the physical matter as a perfect 
fluid with the energy-momentum tensor 


= (p + pWi^U^ + pg,,^, 


(5.1) 


where p = pit), and p = pit) are the pressure, and the energy density of 
the 5 f—matter fluid, respectively. Now, by using the physical matter, and 
the effective energy-momentum tensors together with the FLRW metric Eq. 
fl2.12p in the gf—metric equations Eq. fl2.6p leads us to the t—component 
equation 



as well as the three identical spatial—component equations 


(5.2) 


2a 

a 


fOt\2 k a 2 ~ 

— (-)- :r — 87iGp + mp — 


A3 

Y'- 


(5.3) 


which become the modihed Friedmann equations. By using Eq. (15.21) in Eq. 
(15.31) we can obtain the modihed cosmic acceleration equation as 


d 


AnG 


(3p + p) + 



A3 

Y' 


a 


3 
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(5.4) 
















We observe that the Friedmann, and cosmic acceleration equations are in 
the canonical form, only getting additional contributions from the effective 
fluid pressure, and energy density which are the reflections of the interaction 
Lagrangian term in Eq. fl2.ip which is proportional to the squared graviton 
mass. The matter-fluid equation 

P=- —(p + p), (5.5) 

CL 

is the consequence of the matter energy-momentum conservation law namely, 
~ 0 derived for the FLRW metric. Besides, a similar 

continuity equation 

P = - —(p + p), (5.6) 

CL 

for the effective fluid follows from the substitution of the effective energy- 
momentum tensor Eq. 02.131) into the corresponding Bianchi identity in Eq. 
02.lip upon using the FLRW gf—metric. On the other hand, in this solution 
scheme the /—metric equation becomes 

K [Rl - + p)SVl = (5.7) 

where we have used Eq. 02.191) in Eq. 02.91) , and substituted the result in Eq. 
02.71) . When Eq. 05.6p is satished one does not have to consider the second 
of the Bianchi identities in Eq. 02.lip as it is also automatically satished 
[la m\. Since the two metric sectors are efficiently decoupled from each 
other one can solve the sector equations independently without making 
an assumption on the /—matter. The solution methodology should start 
by hrst choosing which similarity equivalence class representatives in 77i,2,3,4 
links the two metrics. Fixing J in this way determines the equation of state of 
the effective huid via Eq. fl4.ip by substituting the appropriate elementary 
symmetric polynomials of the chosen J. Then, by using the equations of 
state of the effective, and the gf—matter (corresponding to various eras) one 
can solve Eqs. fl5.2p . fl5.5l) . and (15.61) to hnd out the evolution of the scale 
factor, and the state of the effective fluid, and the matter ideal fluid. The 
reader should appreciate that, our solutions are justihed only if one hnds 
also solutions of the /—metric equations Eq. (15.7p . In general, one can 
now use a,p, p (which are completely determined) in Eq. (I4.iup to read the 
associated /—metric which has an implicit dynamical link in Eq. (14.101) to 
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the sector via barely, the metric g, and the elements of J7i,2,3,4 which are 
not only fnnctions of the /5j—parameters of the theory bnt also the effective 
pressnre, and the energy density of the effective ideal flnid whose fnnctional 
forms are solved from the cosmological eqnations of the gf—metric sector. At 
this point we have to remark that, althongh we have previonsly mentioned 
abont a degree of arbitrariness in constrncting / in the set Eq. fl4.12p via 
six arbitrary entries of the matrices Pi these arbitrary fnnctions may also 
be used to £x the form of / now, so that it will satisfy Eq. fl5.7p when one 
chooses the form of /—matter in it. However, this route is not the only one 
to be followed in general. On the contrary, to exemplify a solution outline 
in the /—sector let us consider the solutions A^i, A^ 2 - If they are used in Eq. 
fld.lUp one obtains 



which is in the generalized FLRW form with a lapse function N{t). Here, 
we see that N{t) = (a*)^, and the /—scale factor can be read from b‘^{t) = 
(Ai, 2 )^a^(f). We can also read a* = ai[(3j,p,p], and Ai ,2 = P,P] from 

Eq. 03.191) . Since, as we discussed above from the sector equations a,p,p 
are completely solved N{t), and b(t) in Eq. 05.81) are also determined. Thus, 
in this case the /—metric is hxed as a generalized FLRW one. Let us also take 
the /—matter in the perfect fluid form, and consider = A^{t). With these 
choices, and the substitution of the /—metric from Eq. 05.8p the /—metric 
equation 05.7p becomes 


-m- 


^A^{t)N{t)5^Q5o^ = eSTiG{{pf 

(5.9) 


where 



(5.10) 


We should remind the reader that, in Eq. 05.9p the only unknown functions 
are (t), pf{t),pf{t) as the scale factor b{t) is predetermined. The third 
term on the left hand side in Eq. 05.9p will only contribute a time-dependent 
effective cosmological constant to the 00—component but not to the spatial 
component equations. However, both the 00—, and the spatial-component 
equations will get extra contributions from the lapse function with respect 
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to the FLRW ones. Therefore, from Eq. fl5.9p we will get two modihed 
Friedmann equations which are algebraic (since the scale factor is already 
determined) for the unknown functions A-f(f),There is also a 
hrst-order ordinary differential equation arising from the fluid equation of 
the /—matter perfect fluid, namely, from = 0. Finally, from these 

two algebraic, and one hrst-order ordinary differential equations we can solve 
the unknown functions (t), pf{t),pf{t) to complete the /—sector solution. 


6 Concluding Remarks 

For the massive bigravity theory [9] we constructed the complete solution 
moduli space of the (/, g) couples of metrics which admit a FLRW cosmology 
in the sector via the presence of an effective ideal huid contribution coming 
from the interaction Lagrangian of the mass terms in addition to the matter 
one. We employed the cosmological solution ansatz by choosing the energy- 
momentum tensor of the interaction terms in the metric equations in the 
form of an effective ideal huid one. This choice resulted in a cubic matrix 
equation for the building block matrix of the interaction Lagrangian that 
is composed of the two metrics. By deriving the general solution space of 
this nontrivial matrix equation (whose coefficients are also functions of the 
elementary symmetric polynomials of its solutions) we were able to construct 
and dehne the complete solution space of the (/, g) metric conhgurations 
which enable FLRW cosmologies in the sector that is modihed by an 
ehective ideal huid whose contributions are proportional to the square of the 
graviton mass. Although, we obtained the general solutions of the ansatz 
matrix equation we also discussed that one still has to impose a symmetry 
condition on these solutions to construct a symmetric result for the /—metric. 
Therefore, in spite of the existence of a matrix held degree of freedom in 
constructing /—solutions out of the gf-sector helds one has to render three 
out of nine function components of this arbitrary matrix held to satisfy the 
symmetry condition we mentioned. Furthermore, we also discussed that one 
might also have to hx the remaining degrees of freedom of the /—metric 
in satisfying the /—metric sector held equations in the presence of /—type 
matter. We have shown that, the cosmological solution moduli of bigravity 
that we constructed is composed of similarity equivalence classes which do not 
diher from each other only in their functional form but also in the equations of 
state that they impose on the associated ehective ideal huid they give rise to. 
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Finally, in the last section, we presented the resulting cosmological equations 
of the g, and the /—metrics for which we shortly discussed the solution flow 
chart dictated by the semi-decoupling of the two metric sectors. 

The known exact solutions of bigravity can in general be divided into 
three groups [20]. There is a class of solutions in which both metrics are 
proportional to each other. There exists another class of spherically sym¬ 
metric solutions which has a nondiagonal background metric. There are 
also solutions including both diagonal but not proportional /, and g met¬ 
rics. In this work, we present the complete cosmological background solution 
space of the theory. Massive bigravity as a ghost-free massive gravity theory 
promises to possess cosmological solutions which can exhibit late time self¬ 
acceleration which could compensate the dark energy problem in standard 
cosmology. The background cosmological solutions [121 HSl HH HSl IISI fTT] 
[I81[I91E01[2I1[221[231[211|251E6], and their perturbations and stability issues 
[Ta[22lE3l[30l[3l] arising from the above-listed known solutions have gained 
a considerable interest and they are under extensive inspection recently. It 
has been shown that although there are stability problems and the perturba¬ 
tions of these solutions differ from the ones of GR these problems can still be 
overcome by turning on the /—type matter which is heuristically interpreted 
as dark matter [11122112311301131]. We believe that, apart from its math¬ 
ematical legitimacy of completeness which presents an extensive amount of 
new cosmological solutions of the theory our derivation of the cosmological 
background solution space can also serve for the phenomenology of the the¬ 
ory. We have found that, in the general similarity equivalence class structure 
of the solutions there is a rich variety of functional relations between the 
spatial parts of the two metrics unlike the case in the particular cosmological 
solution which is widely studied in the literature. The behavior of the ratio of 
the g, and /—scale factors of this particular solution (which we believe must 
be related to the Ni, or N 2 solutions we have discussed) causes early time 
instabilities of the perturbations which differ from the GR ones. Therefore, 
we hope that among the variety of complete background solutions we have 
derived there may exist ones which may admit acceptable perturbation be¬ 
havior with respect to GR perturbations. To explicitly construct, and study 
the solutions in this direction one may follow two main routes, one may either 
inspect the solution behavior in the various similarity classes one by one or 
one may attempt to design particular form of cosmological solutions with or 
without /—matter which exhibit a stable nature of perturbations within the 
solution construction methodology we have discussed. However, we should 
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also state that in our generally-constructed solution space, majority of the 
/—metric solutions may fail to exhibit homogeneity, and/or isotropy behav¬ 
ior. On the other hand, one may question the necessity of homogeneity, and 
isotropy in the /—sector since opposite cases may have acceptable results 
from the gf—metric perturbation theory point of view, and in addition they 
may lead to interesting variety of dark matter scenarios. Finally, we point out 
a possible direction in which one can extend the results of the present work 
to study the cosmological solutions within the newly proposed formalism of 
ghost-free effective-metric-matter coupling [32l |33l [M] . 


Appendix 

In the Appendix, we will give a sketch of the proof of the statement we made 
in Section four that any solution of Eq. fl2.19p must belong to the solution 
set fl4.9p . First, let us assume that Eq. fl2.19p has a diagonalizable solution 
such that 

t{Xd)+V = Q. ( 6 . 1 ) 

If we do a similarity transformation which brings to a diagonal Jordan 
form Jd on the above equation then we get 

P-^t{Xd)P + P-^VP = t{Jd) + P~^VP = 0 , ( 6 . 2 ) 

where Jd = P~^X£)P. We can directly observe that P~^PP must be a diag¬ 
onal matrix. Furthermore, since the eigenvalues of P must be invariant under 
a similarity transformation we can conclude that the diagonal matrix P~^PP 
must be in one of the forms; diag(iJ, 0, 0, 0), diag(0, ZJ, 0, 0), diag(0, 0, ZJ, 0), 
or diag(0, 0, 0, D). From this observation we deduce that P must be in one of 
the forms in Eqs. fl4.3|) . or fl4.7l) . This result proves that any diagonalizable 
solution of Eq. fl2.19p must be an element of A4. On the other hand, let 
us consider nondiagonalizable solutions of Eq. fl2.19p with real eigenvalues. 
They also satisfy 

r(XjvD) + ZJ = 0, (6.3) 

which can be brought to a form 

P-^t{Xmd)P + P-^VP = t{Jmd) + P-^VP = 0, (6.4) 
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where = P ^^ndP is one of the nondiagonal Jordan canonical forms 
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where in J 2 , 3 , 4 , 5,6 the distinct diagonal elements 01 , 2 , 3,4 may be equal to e or 
a different value than e. Here, we observe thai 
angular form 

/ Ml U 5 Uq MisX /D 0 0 

U 2 Uq UiQ Mi4 0 0 0 

U 2 , Uj Mil U\Q 0 0 0 

\m4 Ms mi 2 Mie/ Vo 0 0 


p-^VP = 


0\ 

0 

0 

0 / 


t{Jnd) 

is in 

uppertri- 

( 

^2 

^3 

ti\ 


^5 

^6 

h 

ts 


^9 

tio 

til 

tl2 


V^13 

tli 

tl5 

tiej 



[uitiD Mif2-D uit^D uit4^D\ 
U2tlD M2t2-D U2hD U2tiD 
U'itiD M3t2-D u^hD u^tiD 
\uitiD M4f2-D u^t^D UitiD y 


(|6.6p 


must be in uppertriangular form too. Thus, its diagonal elements which are 
its eigenvalues must be {0, 0, 0, D} as under similarity transformations eigen¬ 
values and their algebraic multiplicities are preserved. Therefore, P~^DP 
must be in one of the forms 


(D • • •\ 
0 0 • • 

0 0 0 * 

yo 0 0 0/ 


/o • • •\ 

0 D • • 

0 0 0 * 

\0 0 0 0 / 


/o • • •\ 
0 0 • • 

0 0 T> • 

yo 0 0 0/ 


/o • • .x 
0 0 • • 

0 0 0 • ■ 
yo 0 0 T>/ 
dSZI) 


Let US assume that the nondiagonalizable solution is such that P ^DP be¬ 
comes the second form above. In this case, since M 2 , and ^2 can not be zero 
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we see from Eq. fl 6 . 6 p that U 3 ,U 4 ,ti must vanish. The corresponding Jordan 
form must be Jq because in all the other cases via Eq. flb.jp the diagonal 
entries would lead to two inconsistent equations Ae^ + Be^ + Ce + D = 0, 
and Ae^ + Be^ + Ce = (j^. In this restriction, again from Eq. fl 6 . 6 p we 
see that as t 2 is not zero ui must be zero, also, as U 2 is not zero fa must 
be zero. Besides, Eq. flb.bp denotes that when Jq is used in Eq. fl6.4p 
since M2 7^ 0 ^4 must be zero. Furthermore, P~^P = PP ~^ = 1 yields 
= tio = tii = uq = uio = mi 4 = 0. Therefore, we conclude that in this 
case P = P 2 , and P~^VP = diag(0, ZJ, 0, 0). If we consider the nondiagonal- 
izable solutions which generate the fourth matrix in Eq. fl6.7p for P~^T)P 
we realize that since ^4,^4 are nonzero ^ 1 ,^ 2 , ^3 niust be zero. Also, similar 
to the previous case, since they would result in the inconsistent equations 
Ae^ + Be^ + Ce + D = 0, and Ae^ + Be^ + Ce = 0 the Jordan canonical forms 
J 3 , Jq must be excluded in this case. For this reason, since ^ 0, M 3 must 
be zero as a result of the substitution of the possible Jordan forms J 2 , Ji, J 5 
in Eq. fl6.4p . In addition, since r(J 2 , J 4 , J 5 ) can not have nonzero elements 
at the fourth column except the fourth row, since ^4 7 ^ 0 via Eq. fl 6 . 6 p we see 
that Ml, M 2 must be zero. Upon these substitutions, P~^P = PP~^ = 1 gives 
tg = ti 2 = tiQ = Mg = mi 2 = mi 6 = 0. Hence, we observe that for this case 
P = F 4 , and P~^T)P = diag(0, 0, 0, ZJ). Two straightforward, and similar 
analysis show also that for the hrst, and the third cases in Eq. fl6.7p the trans¬ 
formation matrices must be Pi, and P 3 , also, P~^VP = diag(P, 0, 0, 0), and 
p-^'Dp = diag(0, 0, P, 0), respectively. Besides, the possible Jordan forms 
for these cases are J 3 , 5 , 6 , and J 4 , respectively. We observe also that, Ji is 
not possible for any of the nondiagonalizable solutions with real eigenvalues. 
Therefore, this analysis proves that any nondiagonalizable solution of Eq. 
fl2.19p with real eigenvalues must be contained in A4. Now, let us consider 
the nondiagonalizable soltions of Eq. fl2.19p with complex eigenvalues. By 
applying an appropriate similarity transformation on Eq. fl6.3l) we get 

P-V(Xivd)P + P-^VP = t{Knd) + P-^VP = 0, (6.8) 

where Knd = P~^XndP is one of the nondiagonal Jordan canonical forms 



/Pi 

h 

0 

0 \ 


( R 

z 

0 

0 

\ 

K,= 

-Zi 

0 

Pi 

0 

0 

K2 = 

-I 

p 

0 

0 


0 

P2 

h 

0 

0 

e 

1 



V 0 

0 

-Z2 

P2 / 


VO 

0 

0 

e 

/ 


^We refer the reader to the first footnote in Section three. 
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K3 


( R 
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0 

0 
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( Ol 

0 

0 

o\ 

-I 

R 

0 
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, K^ = 
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R 

I 

0 

0 

0 

03 

0 


0 

-I 

R 

0 

Vo 
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0 

04 

) 


VO 

0 

0 

04/ 


/e 1 0 0\ 

0 e 0 0 

0 0 R I ' 
\0 0 -I R) 


/ Oi 

0 

0 

o\ 

0 

02 

0 

0 

0 

0 

R 

I 

^0 

0 

-I 

r) 


dnsD 


We should state that, when substituted into r all these matrices keep their 
forms with entries changed. For example, 


(R'l 


0 

0 \ 


R'l 

0 

0 

0 

0 

R'2 

R2 

1 0 

0 

-R2 

R'2) 


( 6 . 10 ) 


If we use this in Eq. fl6.8p . and refer to Eq. fl6.6p we see that uRi = ^2^2, and 
U2ti = —Uit2 which give R = t2 = 0 . Also, = ^4^4, and = —u^t^ 
which give ^3 = ^4 = 0 . However, now P becomes a zero-matrix, hence, it 
becomes singular, and can not perform any similarity transformation. There¬ 
fore, this case must be excluded (there exists no nonsingular matrix which 
can bring a nondiagonalizable solution of Eq. fl 2 . 19 p into iFi—form since, 
this results in an inconsistency). Let us consider the case which leads to 
the form 


( a'l 

0 

0 

o\ 

0 

R! 

V 

0 

0 


R' 

0 

^0 

0 

0 

a'J 


( 6 . 11 ) 


Similarly, now, from Eq. fl6.8p we have ^2^2 = ^3^3, and ^2^3 = —^3^2 
which give ^3 = ^2 = 0 . From Eq. fl6.6|) we deduce that, for a consistent 
nontrivial P we can not have ti = 0 , and ^4 = 0 at the same time. Also, 
wi, 2,3,4 can not vanish simultaneously. These facts leave us two cases. Either; 
ti = 0,^4 7^ 0, but Ml, 2, 3 = 0,M4 7^ 0, or R 7^ 6,^4 = 0, but M2,3,4 = 0 ,Mi 7^ 0 
to have consistency when Eq. flb.lip is substituted into Eq. fl6.8p . The hrst 
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case gives P = P 4 , and P~^VP = diag(0, 0, 0, D) (via the preservation of 
the eigenvalnes nnder similarity transformations). Whereas, the second case 
corresponds to P = Pi, and P~^VP = diag(P, 0, 0, 0). Next, let ns consider 
K 2 . Similar to the previons cases above now, for the consistency of Eq. fl 6 . 8 p 
we mnst have = ^2 = 0. Hence, P~^PP mnst be in nppertriangnlar form 
with diagonal elements as its eigenvalnes which mnst be the set {D, 0, 0, 0} 
where D mnst be either at the third, or the fonrth diagonal entry. However, 
for either of these cases Eq. flb.Sp leads us to two inconsistent equations 
Ae^ + Be^ + Ce + D = 0, and Ae^ + Be‘^ + Ce = 0 as we assumed D ^ 0. 
Therefore, this case must be excluded. K 3 on the other hand, leads to the 
conditions Uiti = ^ 2 ^ 2 , and M 2 H = —^ 1^2 that give ^ 1=^2 = 0. Again, Eq. 
flb.bp shows that for a consistent nontrivial P we can not have ^3 = 0, and 
^4 = 0 at the same time, as well as mi, 2 , 3,4 can not vanish all. Thus, either; 
fa = 0 , ^4 7 ^ 0 , but Ml, 2,3 = 0 , M 4 7 ^ 0 , or ts 7 ^ 0 , ^4 = 0 , but mi, 2,4 = 0 , M 3 7 ^ 0 to 
have consistency when rlK^) is used in Eq. flb.Sp . The hrst case gives P = P 4 , 
and P~^VP = diag(0, 0, 0, P), and the second case corresponds to P = P 3 , 
and P~^T>P = diag(0, 0, P, 0). A very similar line of reasoning denotes that 
P 5 is not possible, also, Kq is possible with either; P = Pi, and P“^PP = 
diag(P, 0, 0, 0), or P = P 2 , and P~^VP = diag(0, P, 0, 0). Therefore, we 
conclude that any nondiagonalizable solution with complex eigenvalues of 
Eq. fl2.19p must also be contained in Ad. As a hnal remark, in summary, 
in the Appendix we showed that any diagonalizable or nondiagonalizable 
solution of Eq. fl2.19p must be an element of AA. via proving that their 
Jordan canonical forms must satisfy one of the four equations in Eq. fl2.19p . 
and Eq. fl4.6l) . 
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